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ABSTRACT 

It is often considered that a massive torus with solar mass or so surrounding a stellar- 
mass black hole may be a central engine of a gamma-ray burst. We study the properties 
of such massive accretion tori (or disks) based on the a viscosity model. For surface 
density exceeding about 10 20 g cm -2 , which realizes when about a solar-mass material 
is contained within a disk with a size of ~ 5 x 10 6 cm, we find that (1) luminosity 
of photons is practically zero due to significant photon trapping, (2) neutrino cooling 
dominates over advective cooling, (3) pressure of degenerate electrons dominates over 
pressure of gas and photons, and (4) magnetic field strength exceeds the critical value 
of about 4 x 10 13 G, even if we take 0.1 % of the equi-partition value. The possible 
observable quantum electrodynamical (QED) effects arising from super-critical fields 
are discussed. Most interestingly, photon splitting may occur, producing significant 
number of photons of energy below ~ 511 keV, thereby possibly suppressing pair 
creation. 

Subject headings: accretion, accretion disks — black hole physics — gamma-ray bursts 
— neutrinos 

1. Introduction 



Research on gamma-ray bursts (GRBs) has made rapid progress over the past decade, par- 
ticularly after the discovery of afterglow by Beppo-SAX (see recent reviews by Meszaros, Rees, Sz 
Wijers 1999; Piran 2000; Meszaros 2001). The current common belief is that GRBs are the most 
energetic explosions which ever occurred in the Universe, releasing energy of ~ 10 51 — 10 53 erg only 
in a few to few tens of seconds. It is also widely argued that GRBs result from the conversion of the 
kinetic energy of ultra-relativistic particles created within a fire ball. A relativistic fireball shock 
model was proposed to resolve the compactness problem and has had some success in explaining 
how the GRB and afterglow radiation arise (Rees & Meszaros 1992, 1994; Wijers & Galama 1999). 

However, the central engines of GRBs creating initial hot plasma or extremely energetic par- 
ticles are not well understood yet due mainly to the fact that they are hidden from our view. It 
is usually argued that relativistic phenomena should be somehow involved, since otherwise it is 
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difficult to explain huge fluence and rapid burst-like bust profiles. (Narayan, Paczynski Sz Piran 
1992). Along this line, many interesting possibilities have been proposed so far, such as (1) mergers 
of double neutron star (NS) binaries or NS and BH (black hole) (Paczynski 1986, 1991; Goodman 
1986; Eichler et al. 1989), (2) BH-WD (white dwarf) merger (Fryer et al. 1999); (3) BH-helium 
star merger (Fryer et al. 1999); (4) failed supernova (or collapsars, Woosley 1993; Paczynski 1998; 
MacFadyen & Woosley 1999), and (5) magnetar, rapidly spinning neutron star with extremely large 
magnetic fields (Usov 1992, 1994; Thompson 1994; Spruit 1999). It is of great importance to note 
that almost all the models (except for the magnetar hypothesis) predicts a similar configuration 
as an end result; namely, formation of a few solar mass black hole surrounded by a temporary 
debris torus (or disk) with mass of 0.01 — 1.0 M whose accretion can provide a sudden release of 
gravitational energy (Meszaros, Rees, & Wijers 1999). Then, the key issue will be to understand 
the properties of such a compact and massive disk associated with a huge mass accretion rate. 

Accretion models for gamma-ray bursts were first considered by Narayan, Paczynski & Piran 
(1992), and recently discussed in more details by Popham, Woosley, & Fryer (1999) and Narayan, 
Piran, & Kumar, (2001, hereafter NPK). According to them, neutrino cooling dominates over 
radiative loss at very huge accretion rates (see, however, Chevalier 1996; see also Ruffert &; Janka 
1999 for a numerical simulation of a binary merger including neutrino losses. Therefore, such 
disks (or flows) are called as neutrino-dominated accretion flow (NDAF). In the present study, 
we elucidate the theory of NDAF based on the a viscosity prescription, paying special attention 
to their thermal structure and the implications of expected huge magnetic fields. We present the 
resultant equilibrium solutions and properties of neutrino-cooled disks in section 2. We then discuss 
the effects of huge magnetic fields and implications on neutrino emission in section 3. Final section 
is devoted to conclusions. 



2. Structure of Neutrino-Cooled Disk 
2.1. S-T diagram 

In this study we adopt the following basic equations based on the Newtonian dynamics, for 
simplicity (see, e.g., Kato, Fukue, & Mineshige 1998). We use the cylindrical coordinates, (r,ip,z). 
Then, the mass conservation in the vertical and radial direction is, respectively, represented by 

p = — and M = -2irrv r Y,, (1) 
In 

where p is the matter density, £ is the surface density of the disk, H is the disk half-thickness, M is 
the mass accretion rate, and v r is the radial velocity (negative for inflow). The disk half-thickness 
is, from hydrostatic balance in the vertical direction, given by 



h = Cs /n, 



(2) 
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where the (isothermal) sound velocity c s is defined by 

C s = \1/P, ( 3 ) 

with the pressure p. The angular velocity, £1, is in the Newtonian approximation given by 

n = VgmbhA 3 , (4) 

with G being the gravitational constant and Mbh being the black-hole mass. The pressure is 
composed of three terms: 

P = Prad +Pgas + Pd, (5) 

where p ra d is the radiation pressure, p gas is the gas pressure of the nonrelativistic and non-degenerate 
particles, and p^ is the pressure of the degenerate particles. The radiation pressure is 

2W = ^9*T 4 , (6) 

where g* is statistical degree of freedom in the radiation (=2 for 7, or 11/2 for 7 and e + e~ plasma) 
and a is the radiation constant. The gas pressure is 

Pgas = ^2 riik B T, (7) 

i 

where k-Q is the Boltzmann constant, rij is the number density of particle "i", and the suffix i 
runs the nonrelativistic and non-degenerate particles at the temperature (e.g., e~ and nucleons at 
T< m e c 2 /fee). The pressure of the degenerate particles is written as 

Pd = Pd,nonrcl + Pd,rel- (8) 

Then, 

Pd < nonrci= ^ Air) ~^T' ( ) 

i=nonrcl v a 7 

is for nonrelativistic particles which satisfy the condition of the degeneracy, 

m > giiniikBT^irh 2 ) 3 / 2 , (10) 

for T <C rriiC 2 /kB and m <C rriiC 2 with their mass irn, statistical degree of freedom gi and their 
chemical potential fa, and 

w,.. = Ei(v) 1/3fi2< "' /3 - (11> 

i=rel 

is for relativistic particles which satisfy the condition of the degeneracy, 

ni > gi(k B T/hcf/7r 2 , (12) 
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for T 3> rriiC 2 /k-Q or S> m^c 2 . In this study, for simplicity we assume the complete degeneracy 
whenever the above condition of the degeneracy is satisfied. 

In Fig. 1, we plot the contours of the matter density on the (£, T) plane. From the figure we 
can read off the matter density. It also discriminates the regions where each pressure component 
is dominant and also the regions where electrons and/or nucleons are degenerate. We understand 
that degeneracy pressure is important at high density and low temperature regimes. Such regimes 
inevitably appear in very massive disks, and one can never neglect the contributions by degeneracy 
pressure in GRBs unlike the cases of binary systems or galactic centers. We also notice that one 
should distinguish the regime where only electrons are degenerate from the one where both electrons 
and nucleons are degenerate. 



2.2. chemical potential of electrons 

When electrons are degenerate, there emerges an important consequence on the process of 
neutrino emission, that is, e.g., since the pair creation of electron and positron in the electromagnetic 
thermal bath is suppressed by the charge neutrality (see below), the neutrino emission is also 
suppressed (This is not discussed in NPK). The electron chemical potential fi e is determined by 
the condition of charge neutrality, 

n p = n e --n e +, (13) 

with the number density of electron and positron, 

1 f°° 1 
n e- = I d PP 2 , r^r^ ^ 777^ > ( 14 ) 

n P + 



h 3 TT 2 J e (^/p 2 c 2 +mlc*-^ e )/k B T + 1 ' 



1 f°° 2 1 

Wtt 2 J dpP e (^ P ^cHm^+^)/k B T + \ ' ^ 

In the following sections, it is convenient to introduce a dimensionless chemical potential of electrons 
which is defined by 

r] e = He/k B T. (16) 

Then, the complete degeneracy of electrons corresponds to r\ e 3> 1, and the charge neutrality in 
Eq. (13) becomes 



The chemical potentials of protons and neutrons are defined through the relation, 

roc y 

e (0^ c 2 +m 2 c T_^ )A , BT + 1 



1 f C 1 

= W^Jo fll " ;2 ■ - • (18) 

1 f°° i 
n n = / dpp 2 = . (19) 

H 3 TT 2 Jo e (Vp 2 c 2 +m2 c 4 -Mn)/fc B T + j 
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If the nucleons are not degenerate, we have 

n, = 2(^) 3/2 e-( m ' c2 -" i ^ r , (20) 
for i = p, and n. On the other hand, if nucleons are degenerate, we have 

n * = 3^{—) ' (21) 
for the complete degenerate limit with 

jii = m - niiC 2 . (22) 

In this study, we did not directly solve the exact equation in Eq. (13) to derive the chemical 
potential of electrons. Instead, we performed the following approximated computations in the 
complete-degeneracy regime (rj e » 1). As will be shown later, nucleons are in /3-equilibrium 
at a high temperature T > 5.0 x 10 10 K(Q/0.1) 1 / 5 (r/4r 9 )- 3 / 10 (M B H/3M Q )- 1 / 5 where the neutrino 
emissions are effective and mainly contribute to the cooling process. Then, the chemical equilibrium 
is realized, 

Mp + Me = Mn, (23) 

or 

Mp + Me ~ Q = fin, (24) 

where Q = (m n — m p )c 2 ~ 1.29 MeV. In the regime where only electrons are degenerate, (3- 
equilibrium means 

77 

-^=exp[(Q- /ie )/(fc B r)]. (25) 
n n 

The electron fraction Y e = n v j (n p + n n ) is related with the above ratio by Y e = 1/(1 + (n p /n Tt ) _1 ). 
Combined with the condition of the charge neutrality, i.e., Eq. (17), we may only solve the equation, 

^e x p(-„ e /( tB T)) = ^(Q 3 . (26) 

for ji e » kftT 3> Q. On the other hand, it is relatively easy for us to calculate fi e in the regime 
where both electrons and nucleons are degenerate. From the charge neutrality (Eq. (17)) and 
chemical equilibrium (Eq.(24)), we find that fi e ~ fi n » ftp for k^T <C n e < O(10 2 ) MeV. Then, 
from Eq. (21) we can approximately estimate the chemical potential of electrons as 

Me ^ Mn (27) 

~ 6.628 MeV (p/10 13 g cm" 3 ) 2/3 , (28) 

in the complete-degeneracy limit of both electrons and nucleons. 

We also plot the contours of the chemical potential of electrons in units of MeV (/ie,MeV = 10, 
and 20) in Fig. 1. The chemical potential of electrons play important roles to estimate the emission 
rates of neutrinos which will be discussed in the next subsection. 
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2.3. Heating and cooling rates 

The most important relation in the accretion disk theory is the energy balance between the 
heating and cooling processes. Distinct branches (such as a standard disk or a slim disk, etc.) 
appear because of different heating and/or cooling sources being dominant. According to the 
standard a viscosity model, we express the vertically integrated heating rate (over a half thickness, 
H) as 

Q + = Qvis = l^ 2 , (29) 

where Q+ is denotes the viscous heating rate per unit surface area, the kinetic viscosity v is related 
to the viscosity parameter a by 

v = ^ac s H. (30) 
If we assume the angular momentum conservation, we obtain 

(31) 

From Eqs. (29), (30) and (31), we 

(32) 



i/S 



M 
3tt 



1 




where r m is the radius of the inner edge of the disk (= 3 r g ). 
can relate the heating rate with the mass accretion rate; 



M 




The cooling rate is, on the other hand, summation of three major contributions; 

Q- = Q^A + QZly + Qu, (33) 

where Q~ ad is the radiative cooling rate, Q~ dv is the advective energy transport (Abramowicz 
et al. 1988), and Q~ is the cooling rate due to neutrino loss. Note that instead of including 
advective energy transport in the energy equation NPK adopted the view of CDAF (convection- 
dominated accretion flow; Igumenshchev, Abramowicz, & Narayan 2000, Narayan, Igumenshchev, 
&; Abramowicz, 2000, Quataert & Gruzinov 2000). We, in the present study, retain the classical 
picture based on the (vertically) one-zone treatment (e.g., Kato, Fukue, & Mineshige 1998), since 
it is not yet clear if CDAF provides precise description to the flow structure not only in the low- 
luminosity regimes but also in the hyper-accretion regimes. However, such a distinction is not 
essential here, since we are concerned with the regimes of even higher accretion rates, in which 
neutrino emission is substantial (discussed later). 



The radiative cooling rate is 



(34) 



- 7- 



where a s = ir 2 k B /(60h 3 c 2 ) is the Stefan-Boltzmann constant and the optical depth, r to t, is given 
by 

Ttot = n R pH = (35) 
with being the Rosseland-mean opacity, 

k r = 0.40 + 0.64 x 10 23 (^-3) (|) g- 1 cm 2 . (36) 
The advective cooling rate is given by (Kato, Fukue, & Mineshige 1998), 

Q ^ = ST "^' (37) 
where the radial velocity is v r = —M/(2irrT,) and s denotes entropy per particle, 

s = (s rad + s gas ) / p. (38) 

Here, the entropy density of the radiation is 

Srad = j^<?*T 3 , (39) 



and entropy density of the gas (i.e., nonrelativistic particles) is 

^ , (40) 



/ rp\ 3/2' 



rij V 27r 



where the suffix i runs over nonrelativistic nucleon and electron, gi is the statistical degree of 
freedom of the particle. In the following, we approximated ds/dr as s/r. Note that since the 
entropy of degenerate particles is small in the complete-degeneracy limit, we neglect it here. 

The neutrino cooling rate is composed of four terms; 

Qu = (QNe + <?e+e- + <?brems + <?plasmon)#, (41) 

where q^e is the electron-positron capture rate by a nucleon "N", q e + e - is the electron-positron 
pair annihilation rate, (/brems is the nucleon-nucleon bremsstrahlung rate, and ^piasmon is the rate of 
plasmon decays. Note that NPK considered only the case of the non-degenerate electrons in first 
two terms on the right-hand side of Eq. (41). However, the case of the degenerate electrons in the 
first two terms and the other two terms could be also important. 

The electron-positron capture rate is represented by two terms: 

QNe = Qp+e-^n+Ue + 9n+e+^p+F e ) (42) 
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with 



lp+e -*n+v e 



(1 + Sg A )n p ^ dE e E^El-ml<*{E e - Q) 3 ^^y^T , (43) 
Qn +e+ -, P+ v e = ^3^2(1 + Z9A)n n J ^ dE e E e ^JEI - m'c* (E e + Qf ————{AA) 



2vr 3 /l 3 C 2 
~y2 



where Gp is the Fermi coupling constant (= 2.302 x 10~ 22 cm MeV -1 ), and qa is the axial-vector 
coupling constant of nucleon (~ 1.39) which is normalized by the experimental value of neutron 
lifetime r ra ~ 886.7 s. In the non-degeneracy limit (fi e <C k^T), it is easily estimated by 

*" " 92 x lt,33erg cm " 38 "' (ote) ° (W^*) ■ <45> 

On the other hand, it is a little complicated to estimate the electron capture rate in the electron- 
degeneracy regime. In the only electron-degeneracy regime, we have 

q Ne = 1.1 x 10 31 ?? e 9 erg cm" 3 s' 1 (j^) ' (46) 

for complete degeneracy limit of electrons. Note that in this limit, independently we can also 
obtain the chemical potential of electrons by the balance between the two dominant weak in- 
teraction rates, i.e., q p+e -^ n+l/e = 5.3 x lO 30 ^ 9 . erg cm -3 s^ 1 (T/10 11 K) 9 and q n+e +^ p+Ve = 
2.6 x 10 37 exp(-m e /(/c B T))exp(-7/ e )erg cm" 3 s _1 (p/10 13 g cm" 3 )(T/10 n K) 6 . Compared with the 
chemical potential which was obtained in Eq (26) where we assumed the complete /3-equilibrium, 
they agree with each other just within a few percent, e.g., rj e ~ 3.5 at p = 10 13 g cm -3 and T = 10 
MeV/^B- Therefore, our assumptions of the complete-degeneracy limits would be reasonable. 

When both electrons and nucleons are degenerate, the reactions p+e~ <-> n + v e are suppressed 
by the Fermi blocking of degenerate nucleons in the final states. Then, for example the dE e 
integration in Eq. (43) is limited as \p e , oo], and we approximately obtain the relatively small 
cooling rate, q^ e = 5.0 x 10 32 7yJ erg cm~ 3 s~ - 1 (T/10 11 K) 9 . Compared with Eq. (46), of course, we 
find that this rate is smaller by a factor of l/ry 2 in the complete degeneracy regime, i.e., 7] e S> 1. 

The electron-positron pair annihilation rate through e + + e~ — > v + V is 

q e+e - = 4.8 x 10 33 erg cm'V 1 (j^^j , (47) 

in the non-degeneracy regime, e.g., see Itoh et al. (1989, 1990). When the electrons are degenerate, 
the electron-positron pair annihilation rate is too small to compete the other cooling process, and 
we can neglect it. 

The nucleon-nucleon bremsstrahlung rate through n + n^n + n + v + V is represented by 

*— " 3 ' 4 x 10 " erg cm " 3 s " (i^k) 8 {w^T • (48) 
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in the degeneracy regime of nucleons (Hannestad & Raffelt 1998), and 

= L5 x lo33erg cm " 3 s " (io^k) ' 5 (k^Fs) 2 • 

in the non-degeneracy regime of nucleons (Hannestad & Raffelt 1998; Burrows et al. 2000). 



(49) 



It is also known that the plasmon decay is effective in high densities and high electron degen- 
eracy region (Schinder et al. 1987). The decay rate of the transverse plasmons which are normal 
photons interacting with the electron gas through 7 — > v e + V e is estimated by Ruffert, Janka & 
Schafer (1996) as 

gasmen = 1.5 x 10 32 erg cm" 3 s" 1 (y^) 7 p V>(l + 7 P ) ^ + j-^J , (50) 

where 7 P = 5.565 x lO -2 -^/ (n 2 + 3r/ 2 )/3. In particular, 7 — > u e + V e is a dominant process by a 
factor of ~ 163 compared with the other flavor (— ► v^p^ or v T V T ). 

Note that in NPK only Eqs. (45) and (47) cases were considered in q v . To see what component 
of the cooling rates mainly contributes to the process, in Fig. 2 we plot the most dominant one 
on the (£, T) plane. It is interesting to note that the region where advective energy transport 
dominates covers a rather wide, large X and high T part of the S — T diagram. This is because 
Qadv (oc T 16 ) has a stronger T dependence than Q ra d (oc T 4 ) (shown later). Also notice that 
neutrino cooling is essential only at very high temperatures (T > 10 11 K) and very high matter 
density p > 10 10 - 10 13 g cm~ 3 . The striped region represents the place where heavy elements are 
produced which are expected from the point of view in the equation of state of nuclear matter (Shen 
et al, 1998; Ishizuka, Ohnishi & Sumiyoshi 2002). Here, we adopt the case of the electron fraction 
Y e ~ 0.1 which the numerical simulations of coalescing neutron stars predict, e.g., see Ruffert et al 
(1997). Note that the left side of the boundary of the striped region is just virtual one. Of course, 
the tail realistically continues to the left further. As a result, there exist just a small amount of 
free nucleons. Namely, the neutrino processes which we consider in this study no longer become 
effective there. 

Next, in Fig. 3, we plot the heating and cooling rates as a function of the temperature for 
some representative values of surface density: (a) E = 10 3 ' 5 g cm -2 , (b) S = 10 8 g cm -2 , and 
(c) S = 10 20 g cm -2 . In Fig. 3(a), we find three intersection points. They contain two points in 
which the energy balance is realized due to the radiative or advective cooling processes. As for 
Fig. 3(b), there exists only one intersection point in which the energy balance is realized because 
of the advective cooling processes. It is interesting that in Fig. 3(c), a new branch of the energy 
balance appears at T> m c c 2 //cB, since the neutrino-cooling process becomes effective. By plotting 
only the values at those intersecting points on the (S, T) or (£, M) plane for various values of S, 
we obtain a sequence of the thermal equilibrium solutions; that is the thermal equilibrium curve. 
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2.4. Various time scales 

Here it is necessary for us to check the time scale of the physical processes of the neutrino 
cooling. The dynamical timescale is represented by the accretion time, 



a y GM BH \H 

K OxlO-sm-V^V 72 ^, (52) 

where we assumed that the disk half-thickness is approximately H ~ r/2. On the other hand, the 
time scale of weak interaction of a nucleon off the background electrons and positrons is estimated 
by 

t/3 = {<JN e n e cy l (53) 

~ L0 x 10-4 s (io^k) 5 ' (54) 

where we assume that the cross section is a^e ~ Gp(/cB^) 2 , and the electron number density is 
n e ~ (hsT/hc) 3 . Then, the condition rg < r acc is realized when the following relation is satisfied, 



T > 3.8 x 10 10 K 



/ay/5 / r V 3/1 7MbhV 1/5 



The interaction rate of a neutrino off the background nucleons is roughly given by 

T Nv ~ GK^b^^jvc, (56) 

with the nucleon density tin- Since H ~ r/2, the ratio of the interaction time (1/Tnu) and the 
neutrino crossing time (H/c) is estimated by 

T Nv x H/c * 42.6 ( ( J \ ( ( . (57) 

1 VlO n Ky Vl0 13 g cm- 3 y V4 r g J \3 M Q J V ; 

Since the nucleon to electron-positron ratio is njy/n e ~ O(10)(p/10 13 g cm _3 )/(T/10 11 K) -3 , we 
find that the neutrino interaction is rapid enough to realize /^-equilibrium which we have assumed 
in the previous subsections. When Tn u H/c is greater than unity, we should modify the neutrino 
cooling rate in Eq. (41). To correctly estimate it in that case, we should divide the right-hand 
side of Eq. (41) by Y^ v H/c. However, since the neutrino cooling rate becomes 0(1O 2 ) times larger 
than the heating rate suddenly at intersection points between the cooling and the heating rates, to 
estimate the equilibrium curve, our simple treatments in the previous section are not so wrong and 
the results are not changed. 
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Next, we estimate the vertical diffusion time of neutrinos. It is given by 

Tdiff = — i , (58) 

* Mx , „-^) 2 ( I ^)(^) 2 (^) 2 . (M) 

Compared it with Eq. (51), the condition of Tdiff < face is translated into 

T < 2.3 x 10»K (^)- /2 ( ' 3 ) ~" 2 " 4 (^f) . (60) 

V0.1/ V10 13 gcm- 3 y V 4r <y/ V 3 M ©/ 

Therefore, the neutrino cooling process is effective and realistic in the equilibrium solutions of the 
accretion disk at a high temperature T ~ 10 11 K and a high density p ~ 10 13 g cm -3 which will be 
studied in the next subsection. 



2.5. Thermal Equilibrium Solutions 

According to the procedure discussed in Sec. 2.3, we can find the thermal equilibrium solutions; 
that is, intersecting points of the heating and cooling curves. In Fig. 4 we plot them on the (£, 
T) plane for r = A r g (left panel) and for r = 40 r g (right panel), respectively. In the parameter 
region where neutrino cooling is effective, the dominant processes are the electron capture rate 
(N + e^ N' + u). 

Using Eq. (32), we can also express the mass accretion rate M at the intersection points as a 
function of S. The thermal equilibrium curves on the (£, M) plane are plotted in Fig. 5. Here, 
mass accretion rate is normalized by the critical mass accretion rate, 

M crit = 16L Edd /c 2 , (61) 

where -^Eddis the Eddington luminosity, 

L Ed d = ^GM Bn m v c/a, (62) 

with proton mass m p and the cross section a of the matter and the radiation through the electro- 
magnetic scattering, i.e., -Lgdd/c 2 — 7.3 x 10~ 17 M Q • s _1 (cr/(j T )^ 1 (MBH/-M0) with the Thomson 
cross section ctt — 0.6652 mb. 

There are 5 distinct branches seen in these plots (see Table 1). In the lower-left parts of 
both figures, the equilibrium sequence has an S shape (see Abramowicz et al. 1988), which arises 
because of changes in sources of pressure (gas and radiation pressure) and in sources of cooling 
(radiation and advection). In the upper branch of the S shape, advective energy transport takes 
over radiative cooling. Then, generated photons inside the disk take long time to go out from 
the disk surface so that photons are advected inward and are finally swallowed by a central black 
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hole with accreting matter. Then, disk luminosity is largely reduced. Such trapping effects should 
occur in branches III and IV with even higher M. Neutrino cooling is dominant only in the high 
temperature and density regimes (branch V), in which neutrinos, instead of radiation, can carry 
away the generated energy inside the disk. In branch V of Table 1, we wrote Q~ in the case that 
only electrons are degenerate. On the other hand, in the case that both electrons and nucleons 
are degenerate, Q~ becomes smaller by a factor of 1/r] 2 . Then, the thermal equilibrium solution is 
just realized on the line of the condition of nucleon degeneracy. Namely, the T— £ relation of the 
branch V coincides with that of nucleon degeneracy. In the intermediate branch (IV) still advective 
cooling is dominant, however, note that entropy which is transported by advection is contributed 
by non-degenerate nucleons, since entropy of degenerate particles is little. 

It is interesting to examine the stability of each branch. The criterion for the thermal instability 
whereas that for the viscous (secular) instability is 

(64) 




(see, e.g., Kato, Fukue, & Matsumoto 1998). We understand immediately from Table 1 that only 
branch II (radiation pressure-dominant standard disk branch) is unstable for both modes and all 
other branches, including the NDAF branch, are stable for any modes, as already noted by NPK. 



2.6. Properties of Neutrino-Cooled Disks 

^From the equilibrium solutions found in Sec. 2.5, we derive the following expressions for the 
temperature and density on the equatorial plane and for the corresponding mass-accretion rates as 
functions of radius for a given value of S: 



6/7 / „ \ -9/7 f „ , \ 3/7 



p = 2.3x 10 ia g cm" 3 ( — . ) I — ) ( , ((Hi) 

1 10 20 g cm 2 7 V4r 9 y \S M Q J ' v ' 



and 



M = 1.! X 10 17 Merit W ( 20 S X 1 (-T U f^y 5/ " • («0 
\0.1J Vl0 20 g cm- 2 J \4,r g J \3 M Q J V ; 
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Note that NPK derived similar expressions but in terms of the disk mass, M^k, and the size of 
the initial torus, i? ut, which are roughly related to surface density as 

2M disk 10 20 9 /M disk \ / R out V 2 



7rb Rl ut b 6 V M & J V 5 x 106 cm 

where 6 is a constant of the order of unity, which varies, depending on the geometry of the initial 
torus and we assumed uniform surface density distribution. We thus understand that neutrino 
cooling dominates only if a solar-mass material is contained within a disk with a size of ~ 5 x 10 6 
cm. In other words, if the disk mass is less or if the disk size is larger, neutrino cooling never 
becomes dominant. This constraint is more severe than that obtained by NPK. 

By comparing the left and right panels in Fig. 4, we notice that even at such a high surface 
density as £ ~ 10 20 g cm , the neutrino cooling does not become dominant at a somewhat larger 
radius (i.e., r > 40 r g ). This also supports the claim by NPK that NDAF solution only appears in 
a rather compact region around the center. 

We can get a crude approximation for the strength of magnetic fields, by using the above 
results. The equi-partition field strength, for which magnetic energy is equal to gas energy, is 



-Beq — 



( 8irpk B T \ 1/2 

5/7 / „ X -15/14 , . , \ 5/14 



Therefore, the strength of the magnetic field is very likely to exceed the critical value (e.g. Meszaros 
1992), 

2 3 

B C nt = ^ 4.4 x 10 13 G, (70) 

en 

even if the fields strength is only 0.1 % of the equi-partition value (for which magnetic energy is 
only 10~ 4 % of the particle energy). It is interesting to note that equilibrium B value is about the 
critical field for the proton, 

£crit, P = {m p /m c )B clit ~ 8 x 10 16 G. (71) 

Then various QED (quantum-electrodynamics) effects should manifest themselves (see the next 
section). 



3. Discussion 
3.1. CDAF and slim disk 



The present study is similar to that by NPK, although there are some differences. We newly 
take into account the effects of the electron and nucleon degeneracy in neutrino cooling processes. 
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The resultant disk density and temperature differ in the neutrino-dominated regimes. In the present 
study, we also take the slim-disk picture, instead of the CDAF picture adopted by NPK, at mod- 
erately large £ (and thus M) ranges. As mentioned above, the properties of the hyper-critical 
accretion flow is not yet clear, since full multi-dimensional simulations coupled with radiation hy- 
drodynamics (RHD) and possibly with magnetohydrodynamics (MHD) are finally necessary to 
settle on this issue. Probably convection plays a key role there (Agol et al. 2001), but RHD simu- 
lations so far made in the hyper-critical regimes are based on the approximation of the flux-limited 
diffusion (e.g. Eggum, Coroniti, & Katz 1988; Fujita & Okuda 1998), which may not hold within 
highly turbulent media. Also, magnetic reconnection, which is not easy to simulate in realistic 
situations with extremely high magnetic Reynolds number, seems to play an important role (see 
below) . 

Our conclusion basically agrees with that by NPK in that NDAF appears at small radii only 
for high M disk /r 

out ca ses. NPK concluded that models involving the mergers of BH-NS binaries 
and of NS-NS binaries are favored to explain short GRBs because then the flow is likely to be 
NDAF. However, our condition for the NDAF we found is more stringent that they found and may 
not be generally satisfied. 

NPK also argue that the mergers of BH-WD binaries and BH-helium star binaries are difficult 
to explain GRBs for the following reason. Such mergers will create disks with large dimensions, thus 
being unable to reach the NDAF regimes but leading to the formation of CDAF. Then, convective 
motion mostly transports energy outward, thereby liberating little energy at small radii. It is thus 
unlikely to cause violent explosions like GRBs. However, we would like to point out that the flow 
pattern is rather uncertain in the presence of magnetic fields. 

If there are strong magnetic fields with poloidal (vertical) components, for example, formation 
of outflow (or jet) is unavoidable (e.g. Blandford & Payne 1982; Kudoh, Matsumoto, & Shibata 
1998), and such outflowing matter can produce enormous amount of synchrotron photons. Magnetic 
reconnections, giving rise to intense flare emission, may also occur. That is, mergers of BH and WD 
or helium star cannot be rejected as possible origins of GRBs. More careful discussion is needed 
regarding this issue. 

3.2. Proof of p ra( j > p gas in electron non-degeneracy 

There are some discrepancies in the results between the present study and NPK. Those mainly 
come from a difference in the method. In NPK, for example, they analytically studied the scaling 
law of physical parameters in NDAF in two regimes, i.e., i) the gas pressure-dominated regime, 
and ii) the degeneracy pressure-dominated regime. Here, we show that gas pressure never becomes 
dominant whenever electrons are not degenerate at T> m e c 2 /&b- Namely only when electrons are 
degenerate, either gas pressure or degeneracy pressure can become dominant, i.e., the case i) is not 
realized in electron non-degeneracy. 
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Using Eqs. (6) and (7) with the condition n p ~ n n , we obtain the relation, 




-i 



T 



) 



3 



Prad/Pgas — 



13.02 



P 



(72) 



10 7 g cm 3 



MeV 



On the other hand, the non-degeneracy condition for relativistic electrons, see Eq. (12), is repre- 
sented by 



in the non-degenerate electron regime at T> m e c 2 /kB- 

3.3. Expected QED effects arising due to super-critical magnetic fields 

When the field strength exceeds the critical value (> 0(1O 13 ) G), quantum electrodynamical 
processes become important. Then, as it were, magnetic fields themselves behave as "particles" 
in the strongly magnetized vacuum, thus producing several unique features. For example, strong 
magnetic fields induce large energy splitting of Landau levels, large refractive indices of photons, 
photon splitting effect, and so on (see, e.g., Adler 1971; Shabad 1975; Melrose & Stoneham 1976; 
Chistyakov et al. 1998; Kohri & Yamada 2002a). In particular, the photon splitting effect is 
interesting, when we consider the present model of the hyper-critical accretion disk to a central 
engine of GRBs. That is because if the rate of photon splitting (7 — ► 27) is larger than that of 
e ± pair production (7 + 7 — ► e + + e~), high energy photons with their energy E 1 3>MeV can 
lose kinetic energy rapidly and produce a lot of soft photons (Ej< 0.511MeV) without producing 
copious electron-positron pairs. Then, it may be possible to produce a fire ball as a source of GRB 
without requiring a large Lorentz factor. 

Unfortunately, however, any concrete values of the photon splitting rate near the threshold 
energy of the pair production has not been calculated precisely except for some formal formulations 
(Adler 1971; Shabad 1975; Melrose & Stoneham 1976). The numerical study above the threshold 
is currently underway (Kohri & Yamada 2002b) and needs further attention. 

In addition, radiative processes in a ultra-magnetized plasma are inherently complicated and 
we are still far away from having the basic physics under control. At the enormous densities we are 
dealing with the optical depth is huge, so radiative transfer must be solved. Therefore, the above 
predictions may be premature at present. Only a detailed radiative transfer calculation, which is 
outside the present capabilities of numerical codes, could address this point. 




(73) 



From Eqs. (72) and (73), we see that the following relation is always realized, 



Prad ^ Pgasi 



(74) 
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3.4. Detectability of neutrinos from NDAFs 

As we have shown in the previous sections, we expect that a lot of neutrinos would be emitted 
from NDAF. It is thus of great importance to check whether the signals of the neutrinos are 
detectable or not in GRB events. Moreover, it is also informative to find some correlations between 
neutrino emission and photon emission from GRBs in terms of the observed times and/or directions 
to the source of the neutrinos. 

Such investigation has already been initiated by Nagataki and Kohri (2001), who computed 
the neutrino luminosity and its detectability, properly considering the time evolution of the central 
BH. Their calculations are based not on the merger models but on the collapsar model, and they 
adopt the simple analytical fitting formula to the hyper-critical (slim) disk model taken from Fu- 
jimoto et al. (2001). (Basic numbers are in good agreement with those obtained by the numerical 
evaluations by Popham et al 1999.) Their conclusion is that we will be able to marginally detect 
the neutrinos from collapsars occurring with "TITAND," a next-generation multi-megaton water 
Cherenkov detector (Suzuki 2001). The first phase of TITAND is planed to be with 2 Mt water 
inside. In addition, we will be also able to know the direction to the source of neutrinos in TITAND. 
For the total number of emitted V e , we roughly evaluate 

q Ne VAt 

N v e ~ (75) 

Ev e 

1 n 57 ( QNe \ ( E Ve \ 1 f V \ f At \ 

V10 35 erg cm- 3 s-V V 10MeV / V lO^cnry \ 10" V ' 1 ' 

with their mean energy E Ve to be emitted during one GRB event, where q^e is estimated at 
T ~ lOMeV/fce and p = 10 13 g cm~ 3 , V is the volume of the emitting region, and At is the 
duration. For the event number of F e 's, 

d \~ 2 (N v \ f E Ve \ 2 / Vh 2 o 



R — t ^ 10 {jWc) {wf){wMev) [Jm) (77) 

will be detected by TITAND, where d is the distance from the source and Vh 2 o is the volume of 
water. Hence, we optimistically expect to detect neutrinos from NDAFs, if GRB occurs within a 
distance of 3 Mpc. If detected, neutrino signals and correlations with GRB will be able to make 
clear the unknown explosion mechanism of GRBs. In order to derive the reliable estimations on 
neutrino flux, however, we should perform detailed computations as was done by Nagataki and 
Kohri (2001). This is left as future work. 



4. Conclusion 

We study the properties of hyper-critical accretion flow and for surface density exceeding about 
10 20 g cm -2 , which realizes when about a solar mass material is contained within a disk with a size 
of ~ 5 x 10 6 cm, we find the following unique features: 
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1. Radiation luminosity of such flow is practically zero due to significant photon trapping, al- 
though mass accretion rate enormously exceeds the critical rate, M 3> ^Edd/c 2 - 

2. Neutrino-cooling dominates over advective cooling. Thus the flow can cool via neutrino 
emission. We expect to detect neutrinos from GRBs with the next-generation, multi-megaton 
water Cherenkov detector in future. 

3. Electron degeneracy pressure dominates over gas and radiation pressure, and the degenerate 
electron definitely influences the processes of neutrino emission. This feature is distinct from 
that of the solution found by NPK. 

4. The disks are stable both against thermal and viscous instabilities, even when we take into 
account electron degeneracy. 

5. Magnetic field strength exceeds the critical value, even if we take only 0.1% of the equi- 
partition value. Then, photon splitting may occur, producing significant number of photons of 
energy below ~ 511 keV, thereby possibly suppressing e^ 1 pair creation. However, quantitative 
discussion is left as future work. 

We are grateful to R. Narayan, Shoichi Yamada and S. Nagataki for useful discussions and 
comments. The comments by the anonymous referee are also suggestive and helpful in making the 
final draft. We also thank the Yukawa Institute for Theoretical Physics at Kyoto University for the 
YITP workshops YITP-W-01-13? on "Gamma-Ray Bursts," and YITP-W-01-17 on "Black Holes, 
Gravitational Lens, and Gamma-Ray Bursts," in which this work was initiated and completed. This 
work was supported in part by the Grants-in Aid of the Ministry of Education, Science, Sports, 
and Culture of Japan (13640238, SM). Numerical computation in this work was carried out at the 
Yukawa Institute Computer Facility. 

REFERENCES 

Abramowicz, M.A., Czerny, B., Lasota, J. P., Szuszkiewicz, E. 1988, ApJ, 332, 646 
Adler, S.L. 1971, Ann. Phys, 67, 599 

Agol, E., Krolik, J., Turner, N.J., & Stone, J.M. 2001, ApJ, 558, 543 
Bethe, H.A. 1993, ApJ, 412, 192 

Burrows, A., Young, T., Pinto, P., Eastman, R., & Thompson, T.A. 2000, ApJ, 539, 865 
Chevalier, R. 1996, ApJ, 459, 322 

Chistyakov, M.V., Kuznetsov, A.V., & Mikheev, N.V. 1998, Phys. Lett. B, 434, 67 



-18- 



Eichler, D., Livio, M., Piran, T., k Schramm, D.N. 1989, Nature, 340, 126 
Eggum, G.E., Coroniti, F.V., k Katz, J.I 1988, ApJ, 330, 142 
Fryer, C.L., k Woosley, S. E. 1998, ApJ, 502, L9 

Fryer, C.L., Woosley, S.E., Herant, M., k Daview, M.B. 1999, ApJ, 520, 650 
Fujita, M., k Okuda, T. 1998, PAS J, 50, 639 

Fujimoto, S., Arai, K., Matsuba, R., Hashimoto, M., Koike, O., k Mineshige, S. 2001, PASJ, 53, 
509 

Goodman, J. 1986, ApJ, 308, L47 

Hannestad, S. k Raffelt, G. 1998, ApJ. 507, 339 

Igumenshchev, I.V., Abramowicz, M.A., k Narayan, R. 2000, ApJ, 537, L27 
Ishizuka, C, Ohnishi, A., k Sumiyoshi, K. 2002, nucl-th/0201026 

Itoh, N., Adachi, T., Nakagawa, M., Kohyama, Y., k Munakata, H. 1989, ApJ, 339,354 
Itoh, N., Adachi, T., Nakagawa, M., Kohyama, Y., k Munakata, H. 1990, ApJ, 360, 741 
Kato, S., Fukue, J., k Mineshige, S. 1998, Black-Hole Accretion Disks, Kyoto Univ. Press 
Kohri, K., k Yamada, S. 2002a, Phys. Rev. D65, 043006 

— . 2002b, in preparation 
Kudoh, T., Matsumoto, R., k Shibata, K. 1998, ApJ, 508, 186 
MacFadyen, A. I., k Woosley, S.W. 1999, ApJ, 524, 262 

Meszaros, P. 1992, High Energy Radiation from Magnetized Neutron Stars (U. Chicago Press, 
Chicago), p. 30 

Meszaros, P. 2001, Science, 291, 79 

Meszaros, P., Rees, M.J., k Wijers, R. 1999, New Astron, 4, 313 
Melrose, D.B., k Stoneham, R.J. 1976, Nuovo Cimento, 32, 435 
Nagataki, S., k Kohri, K. 2001, submitted to Prog. Theor. Phys. 
Narayan, R., Igumenshchev, I.V., k Abramowicz, M.A., 2000, ApJ, 539, 798 
Narayan, R., Paczyhski B., Piran T. 1992, ApJ 395, L83 



-19- 



Narayan, R., Piran T., Kumar P. 2001, ApJ 557, 949 (NPK) 

Paczyriski, B. 1986, ApJ, 308, L43 

Paczyriski, B. 1991, Acta Astron., 41, 257 

Paczyriski, B. 1998, ApJ, 494, L45 

Piran, T. 2000, Physics Reports, 333, 529 

Popham, R., Woosley, S.E., & Fryer, C. 1999, ApJ, 528, 356 

Quataert, E., & Gruzinov, A. 2000, ApJ, 539, 809 

Rees, M.J., & Meszaros, P., 1992, MNRAS, 258, p41 

Rees, M.J., & Meszaros, P., 1994, ApJ, 430, L93 

Ruffert, M., & Janka, H.-Th. 1999, A&A, 24, 337 

Ruffert, M., Janka, H.-Th. & Schafer 1996, A&A, 311, 532 

Ruffert, M., Janka, H.-Th., Takahashi, K.& Schafer 1997, A&A, 319, 122 

Shabad, A.E. 1975, Ann. Phys., 90, 1975 

Shen, H., Toki, H., Oyamatsu, K. & Sumiyoshi, K. 1998, Prog. Thor. Phys., 100, 1013. 

Shinder, P.J., Schramm, D.N., Wiita, P.J., Margolis, S.H. & Tubbs, D.L. 1987, ApJ, 313, 531 

Spruit, H.C. 1999, A& A, 341, LI 

Suzuki, Y. 2001, hep-ex/0110005 

Thompson, C. 1994, MNRAS, 270, 480 

Usov, V.V. 1992, Nature, 357, 472 

Usov, V.V. 1994, MNRAS, 267, 1035 

Wijers, R., & Galama, T. 1999, ApJ, 523, 177 



This preprint was prepared with the AAS IATgjX macros v5.0. 



-20- 




logio( S /g' cm " 2 ) 

Fig. 1. — Contours of the matter density on the (S, T) plane for case that the black-hole mass 
is Mbh = 3 Mq and r = 4 r g , where r g is the Schwarzschild radius. The solid lines denote the 
contours of the matter density for log(p/g cm -3 ) = -2, 0, 2, • • -, and 14. The dark (or light) 
shadowed region represents that p gas (p^) dominantly contributes to the total pressure, while the 
white region represents that p Ta d is dominant. The left and right thick dashed lines, respectively, 
denote the loci where electron and nucleon degeneracy begin to take place. Namely, electrons (or 
nucleons) are degenerate because of higher density on the right side of the left (right) dashed line. 
The horizontal dotted line denotes the temperature which corresponds to the electron mass m e . The 
thick dotted line denotes the boundary where electrons become relativistic due to higher energy 
density (p > 2 x 10 6 g cm 3 ) while T < m e c 2 //cB- The dot-dashed lines represent the chemical 
potential of electron in units of MeV, /i e ,MeV = 10, and 20, from the lower left to the upper right. 
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(a) (b) (c) 




5 10 15 20 
logio( 2 /g- cm " 2 ) 

Fig. 2. — Plot of the most dominant component in the cooling rates on the (E, T) plane. The 
adopted parameters are the viscosity parameter, a = 0.1, the black-hole mass, Mbh = 3 M , and 
r = 4 r g . The dark (light) shadowed region and white region, respectively, represent the places 
where the radiative (neutrino) cooling and the advective cooling is the dominant process. On the 
other hand, the region between the two solid lines represents the place where the gas entropy 
s gas mainly contributes to the total entropy. The striped region represents the place where heavy 
elements are produced, and there exist just a small amount of free nucleons. The vertical dotted 
lines, (a), (b), and (c) denote the representative values of E, at which heating and cooling rates are 
plotted in Fig. 3. The dashed and horizontal dotted lines are the same as those in Fig. 2. 
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Fig. 3. — Heating and cooling rates as functions of temperature for a fixed surface density; from 
left to right, (a) £ = 10 35 g cm" 2 , (b) £ = 10 s g cm" 2 , and (c) £ = 10 20 g cm -2 , respectively. 
The other parameters are the same as those in Fig. 2. 
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Fig. 4. — Thermal equilibrium curves (by the thick solid line) on the (£, T) planes at radial 
distances of r = 4 r g (left panel) and r = 40 r g (right panel). The adopted parameters are a = 0.1 
and Mbh = 3 M Q . The vertical dotted lines (a), (b), and (c) denote the three representative values 
of £ adopted in plotting Fig. 3. The thin line indicates the loci where gas pressure equals radiation 
pressure. The dashed line indicates the loci where nucleon degeneracy sets out. We also indicate 
dominant source of cooling by shaded zone as in Fig. 2. In the striped region, there exist just a 
small amount of free nucleons, and the neutrino cooling does not work well. 
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Fig. 5. — Thermal equilibrium curves on the (£, M) plane at r = 4 r g (solid line) , 400 r g (dashed 
line), and 4 x 10 4 r g (long-dashed line). All other parameters are the same as those in Fig. 4. The 
solid circles at the upper ends indicate the places of H = r, above which no physical solution exists. 
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Table 1. T-E and M-£ Relations in Various Regimes 
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'''Here p gas also contributes to the pressure secondly and induces the T-dependence in Q 



